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CAPTIONS 



fig.l : £, £, £ as functions of r. These are called gammal, gamma2, and gammaS, respec- 
tively, on the figures. To make the second and third figures convenient to represent, we 
have added 150 to the values of the second, and 5 x 10*^ to the last, assuring positivity 
(and value greater than unity) in our data range, and taken the logarithm base 10. 
Recursive, approximate quasiperiodic behavior is clearly visible in the variables e, e. 

fig. 2 : The orbit as it enters the region of the attractor, shown in the space of £, £, showing 
the transition from the initial conditions to the strange attractor. In r, the curve 
starts on the e axis to the right of the origin, and ends at the lower left of the origin. 

fig. 3 : This curve continues the evolution within the attractor, starting in the neighborhood 
of the origin, winding up and around to the right and terminating almost at the same 
point (a little above), coming in from the right. 

fig. 4 : This curve is a magnification the motion in the neighborhood of the end point of 
figure 3, starting below the axis and moving to the left. It makes another loop, and 
ends in a smaller pattern. 

fig. 5 : The smaller pattern approached in figure 4 is shown in detail, now displaying a 
somewhat different general pattern. Here, e enters oscillations as e approaches the 
neighborhood of unity. Figures 5a and 5b show enlargements of the innner loop, and 
figure 5c shows the endpoint in more detail. 

fig. 6 : e, e plane in the r interval represented in fig. 2. 

fig7a : £, £ plane in the r interval represented in figure 3. 

fig. 7b : £, e plane in the r interval represented in figure 4, including the details corresponding 
to figure 5. 

fig. 8 : Largest averaged Lyapunov exponent (the values must be divided by the size of the 
interval 10~^) as a function of the accumulated r. We have used interval 0.4 x 10~^. 
Note the sharp jump at approximately 10, 400 steps, indicating the development of a 
very large instability in this neighborhood (the additional terms here are divided by 
10^, the number of time steps. 

fig. 9 : The time series for e{r) on a three dimensional space corresponding to shift by 20 
time steps. 

fig. 10 : The time series for £(r) on a three dimensional space corresponding to shift by 40 
time steps. 

fig. 11 : The time series for eij) on a three dimensional space corresponding to shift by 20 
time steps. 

fig. 12. : In these figures, showing i. vs. e and e vs. £, we have changed the initial conditions, 
and find a similar pattern. 
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Abstract: We review the formulation of the problem of the electromagnetic self-interaction 
of a relativistic charged particle in the framework of the manifestly covariant classical 
mechanics of Stueckeleberg, Horwitz and Piron. The gauge fields of this theory, in gen- 
eral, cause the mass of the particle to change. We show that the non-linear Lorentz force 
equation for the self-interaction resulting from the expansion of the Green's function has 
chaotic solutions. We study the autonomous equation for the off-shell particle mass here, 
for which the effective charged particle mass achieves a macroscopic average value deter- 
mined by what appears to be a strange attractor. 
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1. Introduction 

The advent of powerful computers in the second half of the twentieth century has 
made it a possible and attractive possibility to investigate physical problems that involve 
a high level of intrinsic nonlinearity. The remarkable discovery by Lorenz in 1963 [1] that 
rather simple nonlinear systems are capable of displaying highly complex, unstable, and, 
in many cases, beautiful, systems of orbits, gave rise to the subject of chaos, now under 
intensive study both in specific applications and in its general properties. The instability of 
these orbits is characterized by the fact that small variations of initial conditions generate 
new orbits that diverge from the original orbit exponentially. 

Such studies range from nonlinear dynamical systems, both dissipative and non- 
dissipative, involving deterministic chaos, and hard surface type collisions (for example, 
billliards), to the study of fluctation phenomena in quantum field theory [2]. The re- 
sults have led to a deeper understanding of turbulence in hydrodynamics, the behavior of 
nonlinear electrical circuits, plasmas, biological systems and chemical reactions. It has, 
furthermore, given deep insights into the foundations of statistical mechanics. 

What appears to be a truly striking fact that has emerged from this experience, is that 
the appearance of what has become to be known as chaotic behavior is not just the property 
of some special systems designed for the purpose of achieving some result depending on 
the nonlinearity, such as components of an electric circuit, but seems to occur almost 
universally in our perception of the physical world. The potentialities presented by this 
subject are therefore very extensive, and provide a domain for discovery that appears to 
be virtually unlimited. 

As an example of the occurrence of chaotic behavior in one of the most fundamental 
and elementary systems in nature, we wish to discuss here the results of our study of the 
classical relativistic charged particle (for example, the classical electron without taking 
into account the dynamical effects of its spin, an intrinsically quantum effect). 

There remains, in the area of research on chaotic systems, the question of making a 
correspondence between classical chaotic behavior and the properties of the corresponding 
quantum systems. Much is known on the signatures of chaos for quantum systems, for 
example, the occurrence of Wigner distributions in the spectra for Hamiltonian chaos. The 
study of relativistic phenomena, of the type we shall consider here, where the explicit self- 
interaction problem gives rise to instability, may have a counterpart in the renormalization 
program in quantum field theory. The existence of radiation due to the accelerated motion 
of the particle raises the question of how this radiation, acting back on the particle, affects 
its motion. Rohrlich [3] has described the historical development of this problem, where 
the first steps were taken by Abraham in 1905[4], culminating in the work of Dirac[5], who 
derived the equation for the ideal point electron in the form 

— 

ds"^ ^ ds 

where m is the electron mass, including electromagnetic correction, s is the proper time 
along the trajectory x^^{s) in spacetime, Flf is the covariant form of the electromagnetic 
force tensor, e is the electron charge, and 

_ 2e^ .dV^ _ rf^N 



™-J::^ = ^-''-7r + ^^ (i.i) 
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Here, the indices zv, running over 0,1,2,3, label the spacetime variables that rep- 
resent the action of the Lorentz group; the index raising and lowering Lorentz invariant 
tensor r]^!/ is of the form diag{—l, +1, +1, +1). The expression for was originally found 
by Abraham in 1905 [4], shortly after the discovery of special relativity, and is known as 
the Abraham four-vector of radiation reaction. Dirac's derivation [5] was based on a di- 
rect application of the Green's functions for the Maxwell fields, obtaining the form (1.1), 
which we shall call the Abraham-Lorentz-Dirac equation. In this calculation, Dirac used 
the difference between retarded and advanced Green's functions, so as to eliminate the 
singularity carried by each. Sokolov and Terner [6], for example, give a derivation using 
the retarded Green's function alone, and show how the singular term can be absorbed into 
the mass m. 

The formula (1.1) contains a so-called singular perturbation problem. There is a small 
coefficient multiplying a derivative of higher order than that of the unperturbed problem; 
since the highest order derivative is the most important in the equation, one sees that, 
dividing by e^, any small deviation in the lower order terms results in a large effect on the 
orbit, and there are unstable solutions, often called "runaway solutions". There is a large 
literature [7] on the methods of treating this instability, and there is considerable discussion 
in Rohrlich's book [3] as well (he describes a method for eliminating the unstable solutions 
by studying asymptotic properties). 

The existence of these "runaway solutions", for which the electron undergoes an ex- 
ponential accelaration with no external force beyond a short initial perturbation of the 
free motion, is a difficulty for the point electron picture in the framework of the Maxwell 
theory with the covariant Lorentz force. Rohrlich [8] has discussed the idea that the point 
electron idealization may not be really physical, based on arguments from classical and 
quantum theory, and emphasized that, for the corresponding classical problem, a finite 
size can eliminate this instability. Of course, this argument is valid, but it leaves open 
the question of the consistency of the Maxwell-Lorentz theory which admits the concept 
of point charges, as well as what is the nature of the elementary charges that make up the 
extended distribution. 

It is quite remarkable that Gupta and Padmanabhan[9], using essentially geometri- 
cal arguments (solving the static problem in the frame of the accelerating particle with 
a curved background metric), have shown that the description of the motion of an accel- 
erating charged particle must include the radiation terms of the Abraham-Lorentz-Dirac 
equation. Recognizing that the electron's acceleration precludes the use of a sequence of 
"instantaneous" inertial frames to describe the action of the forces on the electron [10], 
they carry out a Fermi- Walker transformation [11], going to an accelerating frame (as- 
suming constant acceleration) in which the electron is actually inertial, and there solve the 
Coulomb problem in the curved coordinates provided by the Fermi- Walker transformation. 
Transforming back to laboratory coordinates, they find the Abraham-Lorentz-Dirac equa- 
tion without the direct use of the Maxwell Green's functions for the radiation field. This 
result, suggesting the relevance of curvature in the spacetime manifold, such as that gen- 
erated by sources in general relativity, along with other, more elementary manifestations 
of mass renormalization (such as the contribution to the mass due to electromagnetic in- 
teractions and the identification of the Green's function singularity contribution with part 
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of the electron mass), carries an implication that the electron mass may play an important 
dynamical role. 

Stueckelberg, in 1941 [12], proposed a manifestly covariant form of classical and quan- 
tum mechanics in which space and time become dynamical observables. They are therefore 
represented in quantum theory by operators on a Hilbert space on square integrable func- 
tions in space and time. The dynamical development of the state is controlled by an 
invariant parameter r, which one might call the world time, coinciding with the time on 
the (on mass shell) freely falling clocks of general relativity. Stueckelberg [12] started his 
analysis by considering a classical world-line, and argued that under the action of forces, 
the world line would not be straight, and in fact could be curved back in time. He identi- 
fied the branch of the curve running backward in time with the antiparticle, a view taken 
also by Feynman in his perturbative formulation of quantum electrodynamics in 1948 [11]. 
Realizing that such a curve could not be parametrized by t (for some values of t there 
are two values of the space variables), Stueckelberg introduced the parameter r along the 
trajectory. 

This parameter is not necessarily identical to proper time, even for inertial motion for 
which proper time is a meaningful concept. Stueckelberg postulated the existence of an 
invariant "Hamiltonian" K, which would generate Hamilton equations for the canonical 
variables and of the form 

i" = 1^ (1.3) 

and 

where the dot indicates differentiation with respect to r. Taking, for example, the model 
we see that the Hamilton equations imply that 
It then follows that 

- = P (17) 
dt E' ^ ' 

where p^ = E, where we set the velocity of light c = 1; this is the correct definition for the 
velocity of a free relativistic particle. It follows, moreover, that 

. . dx^dXf, p^Pf, 

- - -J^ (1-8) 

With our choice of metric, dx^dx^ = —ds^, and p^p^ = —m?, where m is the classical 
experimentally measured mass of the particle (at a given instant of r). We see from this 
that 

^ = ^ (19) 
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and hence the proper time is not identical to the evolution parameter r. In the case that 
777,^ = M^, it follows that ds = dr, and we say that the particle is "on shell". 

For example, in the case of an external potential V{x), where we write x = x'^, the 
Hamiltonian becomes 

^ = ^ + n^) (1-10) 

so that, since is a constant of the motion, varies from point to point with the 
variations of V{x). It is important to recognize from this discussion that the observable 
particle mass depends on the state of the system (in the quantum theory, the expectation 
value of the operator p^p^ provides the expected value of the mass squared) . 

One may see, alternatively, that phenomenologically the mass of a nucleon, such as 
the neutron, clearly depends on the state of the system. The free neutron is not stable, 
but decays spontaneously into a proton, electron and antineutrino, since it is heavier than 
the proton. However, bound in a nucleus, it may be stable (in the nucleus, the proton may 
decay into neutron, positron and neutrino, since the proton may be sufficiently heavier 
than the neutron). The mass of the bound electron (in interaction with the electromag- 
netic field), as computed in quantum electrodynamics, is different from that of the free 
electron, and the difference contributes to the Lamb shift. This implies that, if one wishes 
to construct a covariant quantum theory, the variables E (energy) and p^ should be inde- 
pendent, and not constrained by the relation E"^ = + where m is a fixed constant. 
This relation implies, moreover, that m? is a dynamical variable. It then follows, quantum 
mechanically, through the Fourier relation between the energy-momentum representation 
of a wave function and the spacetime representation, that the variable t, along with the 
variable x is a dynamical variable. Classically, t and x are recognized as variables of the 
phase space through the Hamilton equations. 

Since, in nature, particles appear with fairly sharp mass values (not necessarily with 
zero spread), we may assume the existence of some mechanism which will drive the parti- 
cle's mass back to its original mass-shell value (after the source responsible for the mass 
change ceases to act) so that the particle's mass shell is defined. We shall not take such 
a mechanism into account explicitly here in developing the dynamical equations. We shall 
assume that if this mechanism is working, it is a relatively smooth function (for example, 
a minimum in free energy which is broad enough for our off-shell driving force to work 
fairly freely) -"^ 

In an application of statistical mechanics to this theory [15], it has been found that a 
high temperature phase transition can be responsible for the restriction of the particle's 
mass (on the average, in equilibrium). In the present work, we shall see that, at least in 



^ A relativistic Lee model has been worked out which describes a physical mass shell 
as a resonance, and therefore a stability point on the spectrum [14] , but at this point it is 
not clear to us how this mechanism works in general. It has been suggested by T. Jordan 
[personal communication] that the definition of the physical mass shell could follow from 
the interaction of the particle with fields (a type of "self-interaction"); this mechanism 
could provide for perhaps more than one mass state for a particle, such as the electron and 
muon and the various types of neutrinos, but no detailed model has been so far studied. 
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the classical theory, the non-linear equations induced by radiation reaction may have a 
similar effect. 

A theory of Stueckelberg type, providing a framework for dynamical evolution of a 
relativistic charged particle, therefore appears to be a natural dynamical generalization of 
the curved space formulation of Gupta and Padmanbhan[9] and the static picture of Dirac. 
We shall study this in the following. 

The Stueckelberg formulation implies the existence of a fifth "electromagnetic" poten- 
tial, through the requirement of gauge invariance, and there is a generalized Lorentz force 
which contains a term that drives the particle off-shell, whereas the terms corresponding to 
the electric and magnetic parts of the usual Maxwell fields do not (for the nonrelativistic 
case, the electric field may change the energy of a charged particle, but not the magnetic 
field; the electromagnetic field tensor in our case is analogous to the magnetic field, and 
the new field strengths, derived from the r dependence of the fields and the additional 
gauge field, are analogous to the electric field, as we shall see). 

In the following, we give the structure of the field equations, and show that the 
standard Maxwell theory is properly contained in this more general framework. Applying 
the Green's functions to the current source provided by the relativistic particle, and the 
generalized Lorentz force, we obtain equations of motion for the relativistic particle which 
is, in general off-shell. As in Dirac's result, these equations are of third order in the 
evolution paramter, and therefore are highly unstable. However, the equations are very 
nonlinear, and give rise to chaotic behavior. 

Our results exhibit what appears to be a strange attractor in the phase space of 
the autonomous equation for the off-mass shell deviation. This attractor may stabilize 
the electron's mass in some neighborhood. We conjecture that it stabilizes the orbits 
macroscopically as well, but a detailed analysis awaits the application of more powerful 
computing facilities and procedures. 

2. Equations of motion 

The Stueckelberg-Schrodinger equation which governs the evolution of a quantum 
state state over the manifold of spacetime was postulated by Stueckelberg[10] to be, for 
the free particle, 

where, on functions of spacetime, is represented by d/dx^ = d^. 

Taking into account full U{1) gauge invariance, corresponding to the requirement 
that the theory maintain its form under the replacement of ip by e'^^°^'ip, the Stueckelberg- 
Schrodinger equation (including a compensation field for the r-derivative of A) is [16] 

(i— + eoa5{x, T))iPr{x) = ^ 

where the gauge fields may depend on r, and cq is a coupling constant which we shall see 
has the dimension £~^. The corresponding classical Hamiltonian then has the form 

^ ^ {p^-eoa'^ix,r))^,-eoa,{x,r)) _ ^^^^^^^^^^ ^2.3) 
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in place of (2.1). Stueckelberg [10] did not take into account this full gauge invariance 

requirement, working in the analog of what is known in the nonrelativistic the 
Hamilton gauge (where the gauge function A is restricted to be independent of time) . The 
equations of motion for the field variables are given (for both the classical and quantum 
theories) by [16] 

Xdaf^{x,T) = eof{x,T), (2.4) 

where cu, /3 = 0, 1, 2, 3, 5, the last corresponding to the r index, and A, of dimension £~^, is 
a factor on the terms f^^fap in the Lagrangian associated with (2.2) (with, in addition, 
degrees of freedom of the fields) required by dimensionality. The field strengths are 

f^P = da^f3 _gP^a^ (2.5) 

and the current satisfies the conservation law [16] 

do,r{x,T) = 0; (2.6) 

integrating over r on {—oo, oo), and assuming that j^{x, r) vanishes at |r| — > oo, one finds 
that 

d^J^ix) = 0, 

where (for some dimensionless rj) [17] 

/CO 
dTi^{x,T). (2.7) 
-oo 

We identify this J^{x) with the Maxwell conserved current. In ref. [18], for example, this 
expression occurs with 

f{x,T) = x^{T)5\x-x{T)), (2.8) 

and r is identified with the proper time of the particle (an identification which can be 
made for the motion of a free particle). The conservation of the integrated current then 
follows from the fact that 

d^^f = x^{T)d^5\x - x{t)) = -^S\x - x{t)), 

a total derivative; we assume that the world line runs to infinity (at least in the time 
dimension) and therefore it's integral vanishes at the end points [10, 18], in accordance 
with the discusssion above. 

As for the Maxwell case, one can write the current formally in five-dimensional form 

= x'^S'^ixir) - x). (2.9) 

For a = 5, the factor x^ is unity, and this component therefore represents the event density 
in spacetime. 
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Integrating the //-components of Eq. (2.4) over r (assuming /^^(x,t) — > for t — > 
±cxd), we obtain the Maxwell equations with the Maxwell charge e = bq/t] and the Maxwell 
fields given by 

poo 

Ai'ix) = A / a^{x,T)dT. (2.10) 



a 



A Hamiltonian of the form (2.3) without r dependence of the fields, and without the 
terms, as written by Stueckelberg [10], can be recovered in the limit of the zero mode of 
the fields (with = 0) in a physical state for which this limit is a good approximation 
i.e., when the Fourier transform of the fields, defined by 

''{x,t)^ j dsa^{x,s)e-''\ (2.11) 

has support only in the neighborhood As of s = 0. The vector potential then takes on the 
form a^(x,r) ~ Asd'^(x,0) = (As/27rA)A'^(x), and we identify e = (As/27rA)eo. The zero 
mode therefore emerges when the inverse correlation length of the field As is sufficiently 
small, and then -q = 27rA/As. We remark that in this limit, the fifth equation obtained from 
(2.4) decouples. The Lorentz force obtained from this Hamiltonian, using the Hamilton 
equations, coincides with the usual Lorentz force, and, as we have seen, the generalized 
Maxwell equation reduce to the usual Maxwell equations. The theory therefore contains 
the usual Maxwell Lorentz theory in the limit of the zero mode; for this reason we have 
called this generalized theory the "pre-Maxwell" theory. 

If such a pre-Maxwell theory really underlies the standard Maxwell theory, then there 
should be some physical mechanism which restricts most observations in the laboratory 
to be close to the zero mode. For example, in a metal there is a frequency, the plasma 
frequency, above which there is no transmission of electromagnetic waves. In this case, if 
the physical universe is imbedded in a medium which does not allow high "frequencies" 
to pass, the pre-Maxwell theory reduces to the Maxwell theory. Some study has been 
carried out, for a quite different purpose (of achieving a form of analog gravity), of the 
properties of the generalized fields in a medium with general dielectric tensor[19]. We shall 
see in the present work that the high level of nonlinearity of this theory in interaction with 
matter may itself generate an effective reduction to Maxwell-Lorentz theory, with the high 
frequency chaotic behavior providing the regularization achieved by models of the type 
discussed by Rohrlich[8]. 

We remark that integration over r does not bring the generalized Lorentz force into 
the form of the standard Lorentz force, since it is nonlinear, and a convolution remains. If 
the resulting convolution is trivial, i.e., in the zero mode, the two theories then coincide. 
Hence, wc expect to see dynamical effects in the generalized theory which are not present 
in the standard Maxwell-Lorentz theory. 

Writing the Hamilton equations 

dx^ dK . dp^ dK . . 

x^=-- = -—,p <mu=-^ = --— 2.12 
dr dr dx^ 

for the Hamiltonian (2.3), we find the generalized Lorentz force 

Mx^' = eor .x" + r 5. (2.13) 
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Multiplying this equation by i^, one obtains 

Mx^x" = eox^ff'^; (2.14) 

this equation therefore does not necessarily lead to the trivial relation between ds and dr 

discussed above in connection with Eq. (1.9). The 5 term has the effect of moving the 
particle off-shell (as, in the nonrelativistic case, the energy is altered by the electric field). 
Let us now define 

e = l+x^'x^ = l- (2.15) 

where ds'^ = dt"^ — rfx^ is the square of the proper time. Since x^ = {p^ — eoa'^)/M, if we 
interpret — eQa^){p^ — eoa^j) = —m^, the gauge invariant particle mass [16], then 

measures the deviation from "mass shell" (on mass shell, ds'^ — dr"^). 

3. Derivation of the differential equations for the spacetime orbit with off-shell 
corrections 

Wc now review the derivation of the radiation reaction formula in the Stueckelberg 
formalism (see also [17]. Calculating the self interaction contribution one must include 
the effects of the force acting upon the particle due to its own field {fseif) in addition to 
the fields generated by other electromagnetic sources (fext)- Therefore, the generalized 
Lorentz force, using Eq.(2.12) takes the form: 

Mx^" = eoX^fea:/ V + CqX'' fself^ v + e^! ext^ 5 + eo/sei/'' 5, (3.1) 

where the dynamical derivatives (dot) are with respect to the universal time r, and the 
fields are evaluated on the event's trajectory. Multiplying Eq.(3.1) by a;^ we get the 
projected equation (2.14) in the form 

M 

—e = eoXf^fext^ 5 + eox^fseif^ 5- (3.2) 

The field generated by the current is given by the pre-Maxwell equations Eq.(2.4), 
and choosing for it the generalized Lorentz gauge ^^a" = 0, we get 

Aa«a"a^(x,T) = {adl - dl + v')a^ = -eo/(x,T), (3.3) 

where cr = ±1 corresponds to the possible choices of metric for the symmetry 0(4, 1) or 
0(3, 2) of the homogeneous field equations. 

The Green's functions for Eq.(3.3) can be constructed from the inverse Fourier trans- 
form 

G{x^ ^) = TTT^ / d^kdn— ^ (3.4) 

^ ' ' (27r)5 J k^k^' + an? ^ ' 



Integrating this expression over all r gives the Green's function for the standard Maxwell 
field. Assuming that the radiation reaction acts causally in r, we shall use here the r- 
retardcd Green's function. In his calculation of the radiation corrections to the Lorentz 
force, Dirac used the difference between advanced and retarded Green's functions in order 
to cancel the singularities that they contain. One can, alternatively use the retarded 
Green's function and "renormalize" the mass in order to eliminate the singularity [6]. In 
this analysis, wc follow the latter procedure. 

The T- retarded Green's function [17] is given by multiplying the principal part of the 
integral Eq.(3.4) by 0{t). Carrying out the integrations (on a complex contour in n] we 
consider the case u = +1 in the following), one finds (this Green's function differs from 
the t- retarded Green's function, constructed on a complex contour in A;°), 



where we have written = x^x^. 

With the help of this Green's function, the solutions of Eq.(3.3) for the self-fields 
(substituting the current from Eq.(2.9)) are 



i^^^^(a;, r) = jj (fx'dT'G{x - x' ,r - r')x^'{r')5^{x' - x{r')) 
= ^ J dr'xi'{r')G{x - x{r'), r - r') 



(3.6) 



and 
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(x, J dT'G{x - x(t'), t - t'). (3.7) 

The Green's function is written as a scalar, acting in the same way on all five com- 
ponents of the source j"; to assure that the resulting field is in Lorentz gauge, however, it 
should be written as a five by five matrix, with the factor 5^ — k^kp/k"^ {k^ = k) included 
in the integrand. Since we compute only the gauge invariant field strengths here, this extra 
term will not influence any of the results. It then follows that the generalized Lorentz force 
for the self-action (the force of the fields generated by the world line on a point x'^{t) of 
the trajectory), along with the effect of external fields, is 

Mx^^^y dT'(x''(T)i;4T')a'^-i"(T)i;'^(T')9.)G(x-x(T'))U=.(r) 

+ i| dr'id^ - x^{r')dr)G{x - a;(r'))U=.(.) ^^'^^ 

We define u = {x^{t) - a;^(T'))(a;''(T) - x^'{T')), so that 

d^ = 2{x^{T)-x^{T'))^. (3.9) 
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Eq.(3.8) then becomes 



Mx'' = 2^ / dT'x''{T)x^{T'){x^'{T) - x^'ir')) 

- x'^iT)x'^iT'){x^{T) - X^iT'))}^G{x - x{t'),T- ^3 .^^^ 

+ dr'(2(a;^(T)-a;^(T0)|^-i:^(T')ajG'(x-a;(T'),T-T')U=.M 

+ eo(/e./.x" + /ext^5)- 

In the self-interaction problem where r — > r', x^{t')—x^{t) — ^ 0, the Green's function 
is very divergent. Therefore one can expand all expressions in r" = r — r' assuming that 
the dominant contribution is from the neighborhood of small t" . The divergent terms are 
later absorbed into the mass and charge definitions leading to renormalization (effective 
mass and charge). Expanding the integrands in Taylor series around the most singular 
point r' = r and keeping the lowest order terms, the variable u reduces to ^ : 

u ^ x^'x^t"'^ - i^xV^ + ix^i V'^ + ^x^x^'t"^. (3.11) 

We now recall the definition of the off-shell deviation, e given in Eq.(2.15), along with 
its derivatives: 

X^x'^ = — 1 + £ 

= 2^ (3.12) 

■■■U .. .. II 1 .. 



Next, we define 



1 . -M^ 1.. 

^^12^^^^ +r 

A = -liT" +wt"\ 



(3.13) 



Using these definitions along with those of Eq.(3.11) and Eq.(3.12) we find 



U + t"2 



T" 



^ In ref. [17], we considered u = x^x^r"^ omitting several possibly significant terms. 
Therefore the corrected equations here contain these additional terms. 
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We then expand the Green's function to leading orders: 

dG ^e{T")hie + A) e{r"),h{e) ef[{e) 



rfiie) £/i(£) : XX 1 ] 



du (27r)V'^ (27r)^ 

dG ^ 29iT")f2{e + A) 2S{T")Me + A) ^ 
dr" (27r)3r"^ {2n)^T"^ 

2e{T") Me) ef^ie) . , ^ ^ 1 -2 w.. xn 1 i^ 



(27r 



|3 



where /' = ^ and 



for £ < 0, 



3tan-^(V=^) ^ ^ 



/2(^) 



3tan-^(v/^) 1 2-s 
(s)^ £2 e2(l-£) 



tan Hy/^) 1 
/3(£) = — -3 + 



(-£)2 



For £ > 0, 



3in|i±v^| 



(£)i ^^(1-^) ^(1-^)' 

ih 

/2(e) = 



IHi^l 1 2-£ 



I £2 £2(1 -£) 



iH^I 1 



£§ £(l-e) 
For either sign of £, when £ ~ 0, 

/2(£)--^-^£-^£'+0(£^), 
/3(e)^^ + ^£+^£^ + 0(£3). 

One sees that the derivatives in Eq.(3.15) have no singularity in £ at £ = 0. 
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From Eq.(3.6) and Eq.(3.7), we have 

fself'' 5{x{t),t) = 

e J {2{x>^{t) - x^{t'))^ - x>^{T')dr}x 

X G{x - x{t'),t - T')\^=^(^r)dr' . 

We then expand x^{t) — x^{t') and x^{t) — x^{t') in power series in r", and write the 
integrals formaUy with infinite hmits. 

Substituting Eq.(3.19) into Eq.(3.2), we obtain (note that a;^ and its derivatives are 
evaluated at the point r, and are not subject to the r" integration), after integrating by 
parts using 5{t") = ^ e{T"), 

M_.^ 2el p gi ^2 £ , ^3 ^ 



2 A(27r)3 



(3.20) 



where we have defined 



gi = fi - f2 - 3/3 , 92 = - h- 2/3 , 93 = g/i - - 2"^^ 
hi = -^fi - 2^2 - /s ) ^2 = - 2/2 - 2"^3 ) ^3 = (/i - /2 - /s) 



(3.21) 



The integrals are divergent at the lower bound r" — imposed by the ^-function; we 
therefore take these integrals to a cut-off /j, > 0. Eq.(3.20) then becomes 

M . _ 2el f gi 92 gs . -v 

2 ^ " A(27r)3 ^ 3^3 -L) 

^i.(e-l)£+^£2 + ^(e-l)«;+ (3.22) 



2//^ 2// // 

+ ^£2(£-l)}^(T")+eoi^/e./5. 
opt, 

Following a similar procedure, we obtain from Eq.(3.8) 

= A&{4«^ - + - f ) + 1^*"^"*" + <^ - 

fi*! .,, 92 93 ■■■u hie... .... h^w . ., h^e'^.,,. (3.23) 

+ f^xi^ - ^x" + —x^ - -^x^ + —x^ + -^x^ + ^— x'^} ^ ^ 

3//-^ 2//^ 2^^ jji S/i 
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Substituting Eq.(3.22) for the coefficients of the terms in the second hne of Eq. (3.23) 
we find: 



• fj, • f V \ ' J 

GqX XiiJ^xt 5 . n a ■ u . f u 

H ; r eo/ea;t vX + e-Qjext 5; 

1 — S 

where 

F(£) = ^(i_e)+^3. (3.25) 

Here, the coefRcients of have been grouped formally into a renormalized (off-shell) 
mass term, defined (as done in the standard radiation reaction problem ) 

M(£, ^) = M + 1^ [ ^'^^~^^ + 92] - [\m -e) + h^js, (3.26) 
where, as we shall see below, 

can be identified with the Maxwell charge by studying the on-shell limit. 

We remark that one can change variables, with the help of (2.15) (here, for simplicity, 
assuming e < 1), to obtain a differential equation in which all derivatives with respect to r 
are replaced by derivatives with respect to the proper time s. The coefficient of the second 
derivative of x^^ with respect to s, and "effective mass" for the proper time equation, is 
then given by 

- ^ [j/i (1 - + + lF(e)VT^ i] 

Note that the renormalized mass depends on £(t); for this quantity to act as a mass, 
£ must be slowly varying on some interval on the orbit of the evolution compared to all 
other motions. The computer analysis we give below indeed shows that there are large 
intervals of almost constant e. In case, as at some points, e may be rapidly varying, one 
may consider the definition Eq.(3.26) as formal; clearly, however, if M{e,s) is large, x'^ 
will be suppressed (e.g., for e close to unity, where Ms{s, e) goes as £/(l — e)^ and M{e, e) 
as s/{l — e)^; note that F{£) goes as (1 — s)~^). 

We now obtain, from Eq.(3.24), 

M{s)x^' = --^^ex^' + F{e)e^x'" + -^x^x"" x^"] 

2 1-^ (3.28) 
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We remark that when one multiphes this equation by x^, it becomes an identity (all 
of the terms except for eofextv^'^ be grouped to be proportional to (^^^^ + ^ii^)', one 
must use Eq.(3.22) to compute the off-shell mass shift e corresponding to the longitudinal 
degree of freedom in the direction of the four velocity of the particle. Eq.(3.28) determines 
the motion orthogonal to the four velocity. Equations (3.22) and (3.28) are the fundamental 
dynamical equations governing the off-shell orbit. 

We remark that as in [17] it can be shown that Eq.(3.28) reduces to the ordinary 
(Abraham-Lorentz-Dirac) radiation reaction formula for small, slowly changing s and that 
that no instability, no radiation, and no acceleration of the electron occurs when it is on 
shell. There is therefore no "runaway solution" for the exact mass shell limit of this theory; 
the unstable Dirac result is approximate for e close to, but not precisely zero. 



4. The £ evolution 

We now derive an equation for the evolution of the off-shell deviation, e, when the 
external field is removed. We then use this equation to prove that a fixed mass-shell is 
consistent only if the particle is not accelerating, and therefore no runaway solution occurs. 
Using the definitions 



Fiie 
F2{e 



3//2 

g2 - 2{e - l)fei 

5'3 + Y^(£- l)/i3 
1 1 

-/i2 + -(£-l)/i4 
1 



(4.1) 



in equation Eq.(3.22), in the absence of external fields, we write 



x^x — 



Fsis) ^2e2 
Differentiating with respect to r we find: 



(4.2) 



X 



c^x'^ + x^x"" = ^{F^e^ + + F2) - F[e - F'^{e)e^ - {2F^{e) + 

F' M 

- F^ie))se - F^s} - -^{Fs + - F^ - F^{s)s' - F^{e)e]e 



(4.3) 



Together with 



(4.4) 
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which is the r derivative of the last equation in Eq.(3.11), one finds from Eq.(4.3), 



Multiplying Eq.(3.28) by (with no external fields) and using Eq.(4.2) and Eq.(4.5), we 
obtain 

(1 + 2§) £ - A{e)e + B{e)e^ + C{e)e - D{e) + E{e)e^ + I{e)ee = 0, (4.6) 

-^3 



where 



2 M 2M{e) AM{e)F, 

Fa ^ ^ e^F{s) 2e^F{e) ' 

„ 2F' ^ M. 2F^ 2 1 .M M{e) 1 

^ = T^i^2-;r:^j -^ + :^ — ^Trl7r:2 +^2) 



F| ^ 2e2^ F3 l-£F3^2e2 e2F(£)l-£' 

4M(£) 1 M 2 2Fi 2 

^ - i^:F^^2i^ + - Tf'^^ - (1 - £)F3 + 

^ 4M(£) Fi 
e2F(£) F3 

(1-£)F3 ^F3^ F3^ 

^^F3+'f3 F3^^ (l-^)i^3- 



(4.7) 



5. Dynamical Behavior of System 

In this section, we present results of a preliminary study of the dynamical behavior of 
the system. The four orbit equations pose a very heavy computational problem, which we 
shall treat in a later publication. We remark here, however, that the quantities (ix/dr and 
dt/dr appear to rise rapidly, but the ratio, corresponding to the observed velocity d'x./dt 
actually falls, indicating that what would be observed is a dissipative effect. 

We shall concentrate here on the evolution of the off-shell mass correction £, and show 
that there is a highly complex dynamical behavior showing strong evidence of chaos. 

It appears that e reaches values close to (but less than) unity, and there exhibits very 
complex behavior, with large fluctuations. The effective mass M{e^ e) therefore appears to 
reach a macroscopically steady value, apparently a little less than unity. In this neighbor- 
hood, it is easy to see from the expression (3.26) that the effective mass itself becomes very 
large, as £/(l — e)^, and therefore the effect of the Lorentz forces in producing acceleration 
becomes small. 

The third order equation of £, Eq.(4.6), can be written as a set of three first order 
equations of the form 

5 = ^(7), (5.1) 
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where 7 = 71,72,78 = £,£,£. We have obtained numerical solutions over a range of 
variation of r. In figure 1, we show, on a logarithmic scale, these three functions; one sees 
rapid fluctuations in 72 and 78, and large variations in 7I. The action of what appears to 
be an attractor occurs on these graphs in the neighborhood of r = 0.002 (time steps were 
taken to be of the order of 10~^). The time that the system remains on this attractor, in 
our calculation about four times this interval, covers about three cycles. The characteristics 
of the orbits in phase space are very complex, reflecting an evolution of the nature of the 
attractor. 

The dynamical signiflcance of these functions is most clearly seen in the phase space 
plots. In flgure 2, we show e = 78 vs. e = 72. This figure shows the approach to the 
apparent attractor inducing motion beginning at the lower convex portion of the unsym- 
metrical orbit which then turns back in a characteristic way to reach the last point visible 
on this graph. We then continue to examine the motion on a larger scale, where we show, 
in figure 8, that this end point opens to a larger and more symmetric pattern. The orbit 
then continues, as shown in figure 4 to a third loop, much smaller, which then reaches a 
structure that is significantly different. The end region of this cycle is shown on larger 
scale in figure 5, where a small cycle terminates in oscillations that appear to become very 
rapid as e approaches unity. The coefficients in the differential equation become very large 
in this neighborhood. 

In figure 6 we see the approach to the attractor in the graph of £ = 72 vs. £ = 7I; 
this curve enters a region of folding, which develops to a loop shown in figure 7a,and in 
fig. 7b, on a larger scale, where we had to terminate the calculation due to limitations of 
our computer. 

We have computed the global Lyapunov exponents according to the standard proce- 
dures; recognizing, however, that there appear to be (at least) two levels of behavior types, 
we have studied the integration on phase space for e somewhat less than unity, since the 
integration diverges for e close to one. For e very close to one, in fact, due to the relation 
ds = (1— £)ciT, this neighborhood contributes very little to the development of the observed 
orbits. 

In this calculation, we have computed the largest Lyapunov exponent (positive) by 
studying the average separation of the orbits associated with nearby initial conditions 
[20]. We show the develpment of the Lyapunov exponent as the average is taken over an 
increasing interval in fig. 8 as data is accumulated up to 18,800 iterations (at least two 
cycles on the attractor). Note that there are intervals of relative stability and very strong 
instability. 

Starting the calulation with different initial conditions, we found very similar dynam- 
ical behavior, confirming the existence of the apparent attractor. 

We have also studied the time series associated with the three variables of the phase 
space, shown in figures 9-11. The autocorrelation functions for e and e fall off quickly 
giving characteristic correlation times; the autocorrelation function for e does not fall off 
faster than linearly, but the time series shows significant structure resembling the phase 
space structure (subtracting out the local average from the function e generates a deviation 
function that does indeed show a fall off of the autocorrelation function) . We have used 
these scales, for 20 time steps (for e and e) and 40 time steps (for e, to contruct a three 
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dimensional delay phase space for the time series for each of these functions. 

In figs. 12a and 12b, we show the behavior of solution of the autonomous equation 
for e vs. e and e vs. e. These results are very similar in form to those given above, 
demonstrating stability of the results. 

6. Discussion 

We have examined the equations of motion generated for the classical relativistic 
charged particle, taking into account dynamically generated variations in the mass of the 
particle, and the existence of a fifth electromagnetic type potential essentially reflecting 
the gauge degree of freedom of the variable mass. We have focussed our attention on the 
autonomous equation for the off-shell mass deviation defined in (2.15) and (2.17) as an 
indication of the dynamical behavior of the system, and shown that there indeed appears to 
be a highly complex strange attractor. It is our conjecture that the formation of dynamical 
attractors of this type in the orbit equations as well will lead to a macroscopically smooth 
behavior of the perturbed relativistic charged particle, possibly associated with an effective 
non-zero size, as discussed, for example, by Rohrlich[8] . Some evidence for such a smooth 
behavior appears in the very large effective mass generated near e = 1, stabilizing the 
effect of the Lorentz forces, as well as our preliminary result that the larger accellerations 
emerging in derivatives with respect to r are strongly damped, in this neighborhood of £, 
in the transformation to motion seen as velocity {d^/dt) in a particular frame. 

Our investigation of the phase space of the autonomous equation used somewhat 
arbitrary initial conditions, and we find that changing these conditions somewhat does not 
affect the general pattern of the results. Moreover, since e depends on the third derivative 
of x^{t), this initial condition should not be taken as arbitrary, but as determined by the 
results for the orbit (initial conditions for the orbit depend only on and its first two 
derivatives). Computing the orbit for a typical choice of initial conditions, and using the 
result to fix e, we find that the corresponding solution of the autonomous equation indeed 
lies in the attractor. The results of our investigation of the autonomous equation should 
therefore be valid for initial conditions that obey the physical constraints imposed by the 
equations of motion as well. 

We have specified all of the dynamical equations here, but a complete investigation 
must await further planned work with more powerful computing facilities and procedures. 
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